Superconductivity due to spin fluctuations originating from multiple Fermi surfaces in 

a double chain superconductor Pr 2 Ba 4 Cu 7 Oi5 ^ 
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The mechanism of superconductivity in Pr2Ba4Cu70is_4 is studied using a quasi-one dimen- 
sional double chain model with appopriate hopping integrals, on-site U, and off-site repulsion V%. 
Applying the fluctuation exchange method to this model and solving the Eliashberg equation, we 
obtain the doping dependence of superconductivity that is consistent with the experiments. The 
superconducting gap has an extended s-wave-like form, which gives a temperature dependence of 
the spin lattice relaxation rate that does not contradict with the experimental results. 
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It goes without saying that unconventional super- 
conductivity (SC) appears on the CuC>2 planes or lad- 
der structures in some cupratesi 2 - Recently, SC in the 
cuprates has been found in another structure, namely, the 
CuO double chain. Pr 2 Ba4Cu 7 0i5_ ( 5 (Pr247/<5), which 
consists of metallic CuO double-chain and semiconduct- 
ing single-chain besides the Mott insulating Cu02 plane, 
shows SC with T™ ax ~ 15 K in a moderate oxygen defect 
concentration range of 6 = 0.2 — 0.6^^, which controls 
the band filling of the double chain block. Recent NQR 
study has revealed that the SC occurs in the double- 
chain, and further observed a "charge freezing" just like 
the one observed in PrBa2Cu40s (Prl24)&£i£ at, and 
only at, 5 = 0, which implies that the double chain is 
near 1/4-filling in Pr247(<5 = 0) (although there remains 
some ambiguity in the band filling estimation) ^ 10 ' 11 

Theoretically, Sano et al. proposed a superconducting 
mechanism using a one dimensional (ID) double chain 
model, 12 which is essentially based on a superconduct- 
ing mechanism proposed generally by Fabrizioi^ and also 
studied by one of the present authors 1 ^ about a decade 
ago. In a purely ID double chain system, the band dis- 
persion can have a double well structure, resulting in four 
Fermi points for appropriate band fillings, which in turn 
results in an opening of the spin gap and dominating 
superconducting correlation. Sano et al. determined the 
tight binding (TB) parameters for a purely ID system us- 
ing the LDA band dispersion between the X and S point 
for the double chain of YBa 2 4 8 (Y124) (See Fig. [2]^. 
From this band structure, they have suggested that the 
appearance of SC upon oxygen reduction is due to the 
increase of the number of Fermi points from two to four 
as the band filling is increased from ~ 1/4-filling. 

Although we believe that this theory is correct in that 
the existence of four Fermi surfaces (FS) is essential in the 
occurrence of SC, there remain several problems, which 
have motivated the present study. First, if we look into 
the band structure in the entire two dimensional Brillouin 
zone (BZ) shown in Figf^l the number of the FS along 
the T-Y line remains to be four down to very low band 
fillings (i.e., the inner FS is 2D) , 16 ' 19 so the variance of 
the FS around 1 / 4- filling is not simply a change of their 
number between two and four as in the pure ID theory. 
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FIG. 1: (a)The effective single band ID double chain 
model, (b) Q1D double chain model. 



Therefore, the reason for the absence of T c before oxygen 
reduction is not clear. Secondly, a previous estimation of 
the spin gap using density matrix renormalization group 
(DMRG) has shown that the spin gap in a double chain 
Hubbard model is, if any, very small at 1/4-filling jii al- 
though the values of the hopping integrals taken there 
do not directly correspond to Pr247. Moreover, a recent 
DMRG analysis shows that the spin gap, if any, is too 
small to be estimated numerically when the ratio between 
the interchain and the intrachain hopping is smaller than 
~ 0.3^ while this ratio is estimated to be around 0.2 in 
Pr247 from first principles calculations. 19 Thirdly, it is 
of interest whether the relatively high T c of 15K can be 
explained within a Hubbard-type model with the above 
mentioned ratio as small as ~ 0.2, because it is obvi- 
ous that SC does not occur when this ratio is too small, 
namely, when the system is essentially a single chain re- 
pulsive Hubbard models To resolve these problems, we 
propose in this paper a spin fluctuation mediated pair- 
ing mechanism, where the origin of the spin fluctuation 
is indeed the presence of multiple FS. We adopt the fluc- 
tuation exchange (FLEX) method to a quasi-one dimen- 
sional (Q1D) extended Hubbard model for the CuO dou- 
ble chain i 21 ' 22 A set of TB parameters are determined 
from the results of LDA calculation for Y124,i^ By 
solving the linearized Eliashberg equation, we obtain a 
finite T c for an extended s-wave SC, whose band filling 
dependence is qualitatively similar to the experimental 
results. 

Note that the original model on the zigzag lattice 
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FIG. 2: (Color online)(a) The TB band dispersion for n = 0.7 
with the determined parameters along with the LDA results 
from ref. flrl (b) The FS (for the noninteracting case) for 
several band fillings. 
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FIG. 3: The obtained T c for singlet pairing. The inset is A as 
functions of temperature. U = 2.0 eV, Vi = are taken. 



(which has two sites in a unit cell and thus results in 
a two band system) can be mapped to a single band 
model having distant hopping integrals as shown Fig. [1] 
Such a mapping is often adopted in the study of zigzag 
chains 11,13 . Note that after this single band mapping, 
the BZ will be unfolded in the kb direction. Hereafter, 
our FLEX results will be shown in the unfolded BZ. The 
Hamiltonian of the extended Hubbard model for the dou- 
ble chain is given by 
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where c. creates an electron of spin a at site i, tij is 
a hopping parameter between site i and j, and U and 
Vij is an electron-electron interaction between such sites. 
The values of Uj are determined by fitting the TB disper- 
sion to the results of the LDA calculations for the dou- 
ble chains of Y124J^ The fitting result is shown in Fig. 
[3J With the obtained values of the hopping parameters, 
we have found that for the noninteracting case the inner 
FS around the T point remains unless n < 0.1 (n=band 
filling=number of electrons/number of sites), although 
a topological change does occur near the X point at 
n ~ 0.3 (Fig. (Hb)), which corresponds to the change of 
the number of Fermi points discussed in ref. \l% In the 
FLEX, the dressed Green's function is obtained by solv- 
ing the Dyson's equation G(k) = [G ( ^" 1 (fc) — E(fc)] -1 self- 
consistently, where Go(fc) is the undressed Green's func- 
tion and S(fc) is the self energy written by the spin and 
charge susceptibilities Xs,c{l) = x(l)[I + VmAtixiq)]' 1 
with k = (k,e n = (2n + l)vrT) and q = (q,u>i = 21ttT). 
x(q) is the irreducible susceptibility which is composed of 
a product of G(fc), and V mi d(q) is the magnetic and den- 
sity coupling vertices described with U and Vij. x(?); 
Xs,c(q) and Vd. m {q) are the matrices indexed by the ini- 
tial (final) relative displacement of particle-hole pair Ar 
(Ar'). The sizes of these matrices increase as the num- 
ber of off-site interactions increases; namely if U and the 
nearest neighbor (n.n.) repulsion V\ exist, the matrix 
size becomes 3x3. The pairing interactions are given by 
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FIG. 4: (Color online) (a) The FS and the nodes of the gap 
function, and (b) the spin susceptibility Xs(l) f° r U = 2.0 eV 
and Vi = eV. 



^VdXcVd](k — k') for triplet pairing. With the obtained 
G(k) and T(k,k'), the linearized Eliashberg equation 
AA(fc) = -^Efe' r ( fc - k')G(k')G(-k')A(k') is solved. 
T c is defined as T c = T(A = 1). In order to assure the con- 
vergence of the calculation, the system size and the num- 
ber of Matsubara frequencies are taken as N — 512 x 16 
sites and 16384, respectively, when only U exists, and 
256 x 16 sites and 16384 are taken when both U and V\ 
are present. As for the band filling range, we have in- 
vestigated in the range of n = 0.3 to near half filling. 
In the following, we show the irreducible susceptibility 
x(<l), spin and charge susceptibilities Xs,c(q) defined as 
the largest eigenvalues of x{q), Xs,c(q)- These along with 
the gap function are plotted at the lowest Matsubara fre- 
quency. 

First, we present the result for [7 = 2.0 eV without the 
off-site interaction. The obtained T c is shown in Fig. [3] as 
functions of n, and the largest eigenvalue A vs. T for var- 
ious band fillings is also shown in the inset. Finite values 
of T c are obtained in a wide range of band filling, and its 
value gradually increases with increasing n. The FS (de- 
fined as £fc— fi+Re(Y,(k)) = 0) along with the nodes of the 
gap function, and the spin susceptibility, are plot- 
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FIG. 5: A vs. T for U = 2.0 eV, Vi = 0.7 eV. (Inset) The 
estimated values of T c obtained by linear extrapolation. 



ted in Fig. Three kinds of nesting (hereafter specified 
by nesting vectors Q 1: Q 2 an d Q3) cause the peaks of Xs- 
The most dominant nesting Q l arises due to the nesting 
between the inner and the outer sets of FS, and thus is a 
direct consequence of the presence of multiple FS, while 
the other two originate from inner-inner or outer-outer 
FS nesting. Q 1 changes from Q 1 — (3?t/4, q a ) at n = 0.5 
to (n/2,q a ) at n = 0.9 and the values of peak itself in- 
crease with increasing the band filling because of the FS 
nesting becomes better. There are no nodes intersecting 
the FS, but the sign of the gap changes its sign between 
the inner and the outer FS, which is a consequence of the 
repulsive pairing interaction mediated by the spin fluctu- 
ations at the nesting vector Q 1 . Consequently, the values 
of A and T c increases with the increase of the peak value 
of Xs at Q x upon increasing the band filling. We will call 
this gap an extended s-wave hereafter. 

The obtained T c of ~ 20 K agrees with the experimen- 
tal result, but the gradual band filling dependence of T c 
fails to explain the "switch on" of the SC with oxygen re- 
duction, i.e., the increase of n in the realistic band filling 
range. At least within the present approach, the change 
of the topology of the inner FS (i.e., change of the num- 
ber of Fermi points along the X-S line) does not strongly 
affect SC. 

This result has motivated us to include the n.n. repul- 
sion V\. Since the system size and the number of Mat- 
subara frequencies are limited, we have calculated A for 
T > 0.004i 2 in order to assure the convergence. Figure 
[5] shows the obtained A as functions of T with V\ = 0.7 
eV. The value of A near n = 0.5 is strongly suppressed 
compared to the result without Vi, and does not seem 
to reach unity even at low temperatures, while for higher 
filling such a tendency of suppression is smaller. This 
result can be traced back to the deformation of the FS 
shown in FigO where the nodes of the gap function and 
the spin and the charge susceptibilities are shown. Al- 
though the positions of Q l and the gap function do not 
differ so much from the results for V\ — 0, the peak value 
of Xs, particularly for n = 0.5, turns out to be smaller. 
As seen in Fig. [7] (a) and (b), where the obtained FS 
and the values of U x(Qi) f° r V% =0 and 0.7 eV are plot- 
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FIG. 6: (Color online) (a) The FS and the nodes of the gap, 
and (b) the spin and the charge susceptibilities Xs and Xc for 
U = 2.0 eV, Vi = 0.7 eV. 



ted, this reduction of Xs is caused by the degradation of 
the nesting condition due to the renormalization of the 
band for V\ ^ 0, which results in a stronger warping 
of the FS especially around 1/4 filling. The reduction 
of Xs is reflected in the reduction of SC especially near 
n ~ 0.5. As for the direct contribution of the charge fluc- 
tuation to the pairing interaction through the ^V^Xc^d 
term, although Xc{Q\) near 1/4 filling is relatively large 
compared to those for other band fillings, the peak value 
itself is much smaller than that of the spin susceptibility, 
so that this contribution should hardly affect SC. Thus 
the main contribution of the charge fluctuation is through 
the self-energy renormalization. 

Although we cannot obtain T c within the investigated 
temperature range, judging the obtained temperature de- 
pendence of A, it does not seem unnatural to linearly ex- 
trapolate log(A) up to A — 1 against log(T) in order to 
estimate the value of T c . In the inset of Fig. [5J the esti- 
mated values of T c are shown. If we assume that the band 
filling of Pr 247/<5 = is n ~ 0.5, this result shows quali- 
tative agreement with the experimental result in that the 
SC appears with moderate electron doping. 

Finally, we discuss the validity of the obtained form 
of the gap function in the light of the experimental re- 
sults. The NQR measurement has observed 1/Ti oc T 2 
(Ti is the spin-lattice relaxation rate) without exhibit- 
ing a coherence peak below T C Z This result seems to 
suggest that the nodes of the gap intersect the FS. How- 
ever, if we assume that a small amount of oxygen de- 
fects controlling the band filling work as scatterers, we 
find that the observed 1/Ti can be accounted for as fol- 
lows. Taking into account the effect of the impurities 
and the defects in the unitarity limit, we have calcu- 
lated 1 jT\ for the gap functions and the FS obtained by 
FLEX*2i As for the temperature dependence of the gap 
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FIG. 7: (Color online) (a) The FS and (b) the filling depen- 
dence of Ux(Qi) for Vi = and 0.7 eV. U is fixed at 2.0 
eV. 



function, A(k,T) = A o 0(fc) tanh(2 v /(T c /T) - 1) with 
A /(A:bT c ) = 4 is assumed, where <j>(k) is proportional to 
the gap function obtained within the FLEX+Eliashberg 
eq. approach. In Fig. [5J the obtained 1/Ti is shown 
for U = 2.0 eV, Vi = 0.7 eV. The values of a/A are 
chosen as 0, 0.2 and 0.4, where a = c/ttNq(0) is the pair 
breaking parameter, c is the impurity concentration and 
iVo(0) is the density of states at the Fermi level in the 
normal state. When a/Ao = 0, 1/Ti decays exponen- 
tially since a full gap opens on the FS. For finite a/Ao, 
1/Ti shows a power-law- like decay without exhibiting a 
coherence peak, which resembles the experimental result. 
Therefore, from these results, we can safely say that the 
observed 1 /Zi does not necessary rule out the possibility 
of the fully gapped extended s-wave pairing. 

To summarize, we have investigated the superconduc- 
tivity of Pr2Ba2Cu40i5-5 using an extended Hubbard 
model on a Q1D double chain with the TB parameters 
determined from the LDA calculation. By applying the 
FLEX to the model and solving the Eliashberg equation, 
in the absence of Vi, we have obtained singlet super- 
conductivity with moderately varying T c in a wide band 
filling range. By introducing the n.n interaction V\, the 
superconductivity is suppressed at a band filling range 
near 1 / 4-filling, which at least qualitatively explains the 
S dependence of superconductivity. NQR relaxation rate 
1/Ti is also calculated taking into account the effect of 
the impurities and the defects. Power-law-like decaying 
1 /T\ resembling the experimental result is obtained with 
a moderate impurity concentration. 

We should comment on the higher 8 regime (6 > 0.5), 



where T c is found to decrease with increasing S and van- 
ishes at S ~ 0.6£ Considering that the nesting vector Q l7 
Q 2 and Q 3 all become close to (w/2,q a ) near half filling 
(see Fig. [Hand Fig. [5]), umklapp processes can be allowed 
to result in a Mott insulating stated which cannot be 
dealt with in the present method. More study on this 
issue will be carried out elsewhere. In a zigzag system 
as in the present case, next n.n. interactions may also 
play some role as suggested in experimental and theoret- 
ical studies.— Unfortunately, the calculation including 
such interactions within the present approach has to be 
restricted to small system sizes. Therefore the study on 
this effect also remains as future study. 
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FIG. 8: The calculated 1/Ti for several values of pair breaking 
parameter for n = 0.7, U = 2.0 eV, Vi = 0.7 eV. 



Present results somewhat resemble that of the recent 
study on Na^OoC^ ■ yH-20 in that the extended s-wave 
superconductivity occurs due to the inner-outer FS nest- 
ing in a disconnected FS system^ Further studies on 
such systems may open interesting new physics. 
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